New generalization of supersymmetric quantum mechanics to arbitrary dimensionality or number of distinguishable particles.

We present here a new approach to generalize supersymmetric quantum mechanics to treat multiparticle and multidimensional systems. We do this by introducing a vector superpotential in an orthogonal hyperspace. In the case of N distinguishable particles in three dimensions this results in a vector superpotential with 3N orthogonal components. The original scalar Schrödinger operator can be factored using a 3N-component gradient operator and introducing vector "charge" operators: Q(1) and Q(1)(dagger). Using these operators, we can write the original (scalar) Hamiltonian as H(1) = Q(1)(dagger) x Q(1) + E(0)((1)), where E(0)((1)) is the ground-state energy. The second sector Hamiltonian is a tensor given by H(2) = Q(1)Q(1)(dagger) + E(0)((1)) and is isospectral with H(1). The vector ground state of sector 2, psi(0)((2)), canbe used with the charge operator Q(1)(dagger) to obtain the excited-state wave function of the first sector. In addition, we show that H(2) can also be factored in terms of a sector 2 vector superpotential with components W(2j) = -(partial partial differential ln psi(0j)((2)))/partial partial differentialx(j). Here psi(0j)((2)) is the jth component of psi(0)((2)). Then one obtains charge operators Q(2) and Q(2)(dagger) so that the second sector Hamiltonian can be written as H(2) = Q(2)(dagger)Q(2) + E(0)((2)). This allows us to define a third sector Hamiltonian which is a scalar, H(3) = Q(2) x Q(2)(dagger) + E(0)((2)). This prescription continues with the sector Hamiltonians alternating between scalar and tensor forms, both of which can be treated by the variational method to obtain approximate solutions to both scalar and tensor sectors. We demonstrate the approach with examples of a pair of separable 1D harmonic oscillators and the example of a nonseparable 2D anharmonic oscillator (or equivalently a pair of coupled 1D oscillators). We consider both degenerate and nondegenerate cases. We also present a generalization to arbitrary curvilinear coordinate systems in the Appendix.